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ABSTRACT 

The elliptical galaxy NGC 1550 at a redshift of z = 0.01239, identified with 
an extended X-ray source RX J0419-I-0225, was observed with XMM-Newton for 
31 ks. From the X-ray data and archival near infra-red data of Two Micron All 
Sky survay, we derive the profiles of components constituting the NGC 1550 sys- 
tem; the gas mass, total mass, metal mass, and galaxy luminosity. The metals 
(oxygen, silicon, and iron) are extended to ~ 200 kpc from the center, wherein 
~ 70% of the i^-band luminosity is carried by NGC 1550 itself. As first revealed 
with ASCA, the data reconfirms the presence of a dark halo, of which the mass 
(1.6 X 10^^ Mq) is typical of a galaxy group rather than of a single galaxy. Within 
210 kpc, the K-hand mass-to-light ratio reaches 75 Mq/ Lq, which is comparable 
to those of clusters of galaxies. The iron-mass-to-light ratio profile (silicon- and 
oxygen-mass-to-light ratio profiles as well) exhibits about two orders of magni- 
tude decrease toward the center. Further studies comparing mass densities of 
metals with those of the other cluster components reveal that the iron (as well 
as silicon) in the ICM traces very well the total gravitating mass, whereas the 
stellar component is significantly more concentrated to within several tens kpc 
of the NGC 1550 nucleus. Thus, in the central region, the amount of metals is 
significantly depleted for the luminous galaxy light. Among a few possible ex- 
planations of this effect, the most likely scenario is that galaxies in this system 
were initially much more extended than today, and gradually fell to the center 
and merged into NGC 1550. 

Subject headings: galaxies: clusters: individual (NGC 1550) — galaxies: elliptical 
and lenticular, cD — X-rays: galaxies: clusters 
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1. Introduction 

Early-type galaxies, despite their optical homogeneity (Kodama & Matsushita 2000), 
exhibit two distinct subclasses with respect to their X-ray properties (Matsushita 2001; 
Vikhlinin et al. 1999); X-ray extended galaxies, and X-ray compact ones. Compared to 
objects in the latter class, those in the former class are more luminous and extended in X- 
rays, and are surrounded by group-scale potential halos (Matsushita et al. 1998; Kawaharada 
et al. 2003). Some of such X-ray extended ones are cD or XD (X-ray dominant) galaxies 
in groups or clusters, while others appear optically isolated, which we call isolated X-ray 
extended galaxies (IXEGs). 

IXEGs are also called fossil groups (Jones et al. 2003), or X-ray overluminous elliptical 
galaxies (Vikhlinin et al. 1999). The original definition of a fossil group (Jones et al. 2003) 
is that its X-ray bolometric luminosity is higher than 4.8 x 10^^ erg s~^ (for a Hubble 
constant of 72 km s~^ Mpc~^), and that the two brightest member galaxies within half 
the virial radius have an i?-band magnitude difference (Ami2) larger than 2. Jones et al. 
(2003) thus identified five fossil groups, and found that they have significantly higher X-ray 
luminosities than normal groups of similar total optical luminosities. These objects were 
thought to represent 8-20 % of all systems with the X-ray luminosity of > 4.8 x 10^^ erg 
s~^. Later, more than 30 objects have been identified as fossil group candidates (Santos 
et al. 2007) by combining the Sloan Digital Sky Survey (SDSS) data with those from the 
Rosat All Sky Survey (RASS). From a study of seven fossil groups, Khosroshahi et al. (2007) 
found that they show systematically higher X-ray luminosities and temperatures than normal 
galaxy groups of similar optical luminosities, and higher mass concentrations. These authors 
attributed the more enhanced X-ray emission from these fossil groups to their cuspy potential 
shape. 

As suggested by cosmological simulations (D'Onghia et al. 2005; Dariush et al. 2007), 
a likely formation scenario of an IXEG (or a fossil group) is that it was an ordinary galaxy 
group at the beginning, and its member galaxies have merged into a central galaxy which 
is now observed as an IXEG (e.g. Jones et al. 2003). Indeed, a rich population of dwarf 
galaxies, found around the IXEG NGC 1132 (Mulchaey & Zabludoff 1999), provides good 
evidende, because such dwarf members are expected to have very long merging time scales, 
whereas the brightest group members are thought to have merged in a few tenths of the 
Hubble time if the galaxy density was initially high enough (e.g. Barnes 1989; Zabludoff & 
Mulchaey 1998). Nevertheless, we may not be able to rule out an alternative possibility that 
IXEGs have been one-galaxy systems since their formation, and their group-like properties 
are inherent to their birth process. In the present paper, we address this issue from X-ray 
observations of one such object, NGC 1550. 
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Among IXEGs ever found, the SO galaxy NGC 1550 is one of the nearest (at a redshift 
oi z = 0.01239) and the X-ray brightest. Historically, an extended (~ 14' in radius) X-ray 
emission named RX J0419-I-0225 was first discovered by the RASS from a position centered 
on this galaxy. Although the NGC 1550 system does not completely meet the condition of 
Ami2 > 2 for fossil groups given by Jones et al. (2003), it has only a few other faint galaxies 
within the bright extended X-ray emitting region, and hence an ideal target to explore the 
nature of IXEGs. In 1998, we observed this object with ASCA for 48 ks; assuming the 
Hubble constant as Hq = 72 km s~^ Mpc~^, we found that the temperature (1.4 keV), X-ray 
luminosity (5.1x10"^^ erg s~^), gas mass (3.4x10^^ ^o), and the total mass (LGxlO"*^^ Mq) of 
this object are typical for those of galaxy groups, rather than of single galaxies (Kawaharada 
et al. 2003; hereafter Paper I). In addition, its mass-to-light [B band) ratio turned out to be 
as high as ^ 380 Mq/Lq, comparable to those of galaxy clusers. Thus, NGC 1550, associated 
with the X-ray source RX J0419-I-0225, can be regarded as a typical IXEG, surrounded by 
a group-scale gravitational potential. 

With its high surface brightness in X-rays and its large X-ray angular extent, NGC 1550 
is expected to provide a new X-ray test of the galaxy merger scenario. It is widely agreed 
that metals in the X-ray emitting plasmas (Intra-Cluster Medium, or ICM) in galaxy groups 
and clusters have been supplied by member galaxies, either at an early stage of the system 
formation, and/or over the Hubble time (e.g Matteucci & Vettolani 1988). Furthermore, 
the metals ejected to the ICM are not expected to diffuse significantly from the place where 
they were first deposited. For example, an iron ion ejected into the ICM from a galaxy is 
slowed down by collisions with surrounding protons and helium ions (Rephaeli 1978). Since 
the diffusion constant of iron in the ICM is ~ 10^^ cm~^ s (Ezawa et al. 1997; 2001), iron 
ions can diffuse at most ~ 10 kpc over the Hubble time, assuming no bulk flows. Then, if 
NGC 1550 is a one-galaxy system from the birth, its metals should concentrate within the 
region of the central galaxy, ~ 30 kpc in radius. Otherwise, the metals should extend over 
the group-scale potential well. Since the ASCA data were too limited in statistics to reliably 
answer this inquiry, we resorted to XMM-Newton observations. Our proposal based on this 
idea was accepted in the XMM-Newton AO-2 cycle, and the observation was carried out 
successfully in 2003. The prime goal of the present study is to compare spatial distributions 
of metals in the ICM and the stellar mass. We utilize the archival near infra-red data of 
Two Micron All Sky Survey (2MASS) for the information of galaxy light. 

We report the XMM-Newton observation of NGC 1550 in §2, followed by the analysis 
and results of the EPIC data in §3. In §4, we briefly evaluate the near infra-red data of 
NGC 1550 and surrounding galaxies. The X-ray and near infra-red data are compared in 
§5, focusing in particular whether the metals in the ICM are spatially more extended than 
the stellar light. In §6, we discuss the nature of the NGC 1550 system based on these 
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results. Quantities quoted from the literature are rescaled to Hq = 72 km Mpc~^. The 
definition of solar abundance is from Anders and Grevesse (1988). Errors are given at the 
90 % confidence level. 



2. Observation 

2.1. Data preparation 

NGC 1550 was observed with XMM-Newton on 2003 February 22 for 31 ks, with the 
Prime Full Window mode and medium filter for the EPIC, and the High Event Rate with 
SES mode for the RGS. To create event files, we processed the Observation Data File of the 
PN and MOS data using "epchain" and "emchain" tasks of the standard analysis package 
(SAS) version 6.5.0, respectively. The created event files were further filtered with respect 
to the event quality fiag and pixel patterns. We chose events which have good quality (not 
collected at the edge of CCDs or bad pixels) and single pixel pattern for PN, and good 
quality and less than quadruple pattern for MOS taking into account the ~ 4 times smaller 
MOS pixel size (l."l squared). In order to reject background fiares, we made 10-12 keV 
count-rate histograms of each detector with 100 sec integration, and discarded those time 
bins which are deviated by more than 2a from the mean rate. After this event selection, the 
effective exposure became 19.6 ks (PN), 21.1 ks (MOSl), and 21.4 ks (M0S2). Point sources 
were excluded by "ewavlet" task of the SAS with the detection threshold set at 7a. Then, 
we made response files for each CCD through "rmfgen" and "arfgen" tasks. 

2.2. Background estimate 

For the present purpose, an accurate background subtraction is of essential importance: 
in order to measure the ICM metallicity out to large radii where the X-ray surface brightness 
becomes low, we need to accurately determine not only atomic lines, but also contnuum. As 
template background files of the three EPIC detectors, we used those data taken with the 
medium filter, publicly released by the University of Birmingham group0 (Read and Ponman 
2003). These background files are superpositions of multiple observations, and hence have 
higher statistics and smaller systematical errors than other background datasets currently 
available. These background events were rearranged to match the observational condition 
of NGC 1550, as originally developed with the ASCA CIS (Ishisaki 1995) and also applied 
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to the XMM-Newton EPIC data by Kaastra et al. (2004). That is, we sorted the entire 
background data into subsets according to 100-sec averaged 10 — 12 keV count rates, and 
created background spectra separately for individual subsets. We then combined them into 
a single background spectrum (one for each detector), using appropriate weights which are 
specified by the 10 — 12 keV count rate histograms (averaged over 100 sec) of the actual 
NGC 1550 observation. For example. Figure [1] is the MOSl 10 — 12 keV count (per 100 
sec) distribution during the NGC 1550 observation, compared with that of the template 
background. Thus, the two distributions are reasonably matched with each other. The 
background spectrum to be utilized in the present analysis, B{E), is then synthesized as 

B{E) = ^^ME) 

where E is the energy, Tj is the total residence time when the 10 — 12 keV count is in the 
j-th interval, and Bj [E) is the background template for the same j-th interval. 

Katayama et al. (2004) studied the properties of the PN background and found that 
2 — 7 keV background rate correlates tightly with that in the 10 — 12 keV band. Its correlation 
coefficient is 0.929, and the scatter around the correlation is less than 8%. Kaastra et al. 
(2004) showed that the MOS 0.2 — 10 keV background flux scales proportionally to the 
10 — 12 keV background rate, and the spectral shape does not depend on the rate within 
an uncertaintly of 10%. Therefore, the above method of background synthesis is considered 
to reproduce the background with an accuracy better than 10%. Figure [2] compares the 
0.5 — 9.0 keV raw EPIC spectra of NGC 1550 with the estimated background spectra. In 
the higher energy range above ~ 6 keV, the two spectra agree within 10% (lo"). 

Since the Galactic absorption toward NGC 1550 is rather high as 1.15 x 10^^ cm~^ 
(Dickey & Lockman 1990), the cosmic X-ray background (CXB) which dominates the low 
energy band might be overestimated when the summed blank sky background is used. How- 
ever, at 0.6 keV, the count rate ratios among the signal, particle background, and CXB is 
estimated to be 1 : 0.2 : 0.06. If the Galactic absorption is 3 x 10^° cm~^, this ratio becomes 
1 : 0.2 : 0.09. Therefore, the difference of the CXB level caused by using the averaged blank 
sky background is at most ~ 3% and ~ 15% of the signal and particle background at 0.6 
keV, respectively. In order to deal with this uncertainty of cosmic X-ray background, and of 
the particle background as well, we put systematic errors to the spectra, following Kaastra 
et al. (2004). The systematic errors in the source spectrum were set at 5% in 0.5 — 2.0 keV, 
and 10% above 2.0 keV. To the background spectrum, we assigned errors of 25% in 0.5 — 0.7 
keV, 15% between 0.7 — 2.0 keV, and 10% above 10 keV. We added these systematic errors 
in quadrature to the corresponding statistical errors. 
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3. X-ray Data Analysis and Results 

3.1. Radial surface brightness 

As shown in Figure [3] (left), the X-ray emission of NGC 1550 is spatially round, and 
extends at least up to a 2- dimensional radius of 6 ~ 14' (210 kpc). The X-ray emission 
centroid coincides with the position of the NGC 1550 nucleus within 10" (~ 3 kpc). The 
signal-detected region is consistent with that obtained by ASCA (Paper I). 

Figure H] shows the azimuthally-averaged 0.5 — 4.5 keV radial surface brightness pro- 
file S(6), taken with MOSl and presented after subtracting the blank-sky background as 
prepared above. We fitted the surface brightness profile by a /9 model convolved with point- 
spread function (PSF). The position- and energy-dependence of the PSF was quoted from the 
XMM-Newton calibration note EPIC-MCT-TN_011 (Ghizzardi 2001). The single-/? model, 
however, did not reproduce the radial profile well, with x^/^ = 1446/599 (/3 = 0.45, and the 
core radius of = 5 kpc). Therefore, we added another j3 componet to the model as 

^ -3/3i+l/2 f , , 9^ -3/32+1/2 



where Si and S2 are normalization constants, while ri and r2 are core radii. This "double-/?" 
model has been successful, yielding x^/ ^ = 677/598. The inner (3 component is characterized 
by /3i = 0.60 lb 0.05 and a core radius of ri = 3 ± 1 kpc, while those of the outer component 
are /?2 = 0.50 ± 0.05 and r2 = 20 ± 3 kpc. These two components are shown separately in 
Figure HI 

In Paper I, a single-/? fit to the ASCA data was successful, and yielded /? = 0.47 ± 0.01 
and = 16 ± 1 kpc. The value of (3 derived with ASCA is consistent with that of the 
outer /^-component (/?2) obtained here, while the core radius of ASCA falls between ri and 
r2 of XMM-Newton. We interpret that the angular resolution of ASCA was not sufficient to 
resolve the inner-/? component, which is seen within ~ 15 kpc or ~ 1'. 

According to Sun et al. (2003), the surface brightness profile of NGC 1550 obtained 
with Chandra cannot be reproduced by a single-/? model either, and requires a double-/? 
modeling. The derived outer component has /? = 0.48 and rc ~ 26 kpc, while the inner one 
is not described in their paper. Thus, the double-/? property, found with XMM-Newton, is 
also supported by the Chandra data. Sun et al. (2003) also found an excess emission within 
1 kpc even above the double-/? fit; this is also suggested in our result, although the angular 
scale of 1 kpc, or ~ 4", is comparable to the PSF of XMM-Newton. 
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3.2. Spectral analysis 

In order to investigate global spectral properties of NGC 1550, we collected photons over 
a circular region of 14' in radius from the galaxy nucleus. The background was subtracted 
in the same way as described in § 2, applying the same data integration region to the 
background dataset as the on-source data accumulation. The spectra obtained with the 
three EPIC detectors were fitted with an absorbed vMEKAL model (Kaastra & Mewe 1993) 
or an absorbed vAPEC model (Smith et al. 2001), using XSPEC 11.3.2. To secure a stable 
fit convergence to physically reasonable parameter values, we fixed abundances of He, C and 
N to one solar, constrained those of other metals from zero to 3 solar, and fixed the column 
density to the Galactic value (1.15 x 10^^ cm~^). The employed energy range is 0.5 — 5.0 keV 
for MOS, and 0.7 — 5.0 keV for PN. The 0.5 — 0.7 keV energy region of PN was excluded, 
because of systematic residuals below Fe-L lines often found in PN spectra of bright clusters 
(Tamura et al. 2004). 

We have found that a reliable estimation of oxygen abundance requires a particular 
caution. This is because the oxygen lines of NGC 1550 are weak, and the model con- 
tinuum fitted over the 0.5 — 5.0 keV range (or 0.7 — 5.0 keV for PN) often becomes slightly 
(~ 5%) higher (or lower) around the oxygen lines than the true continuum, probably due to 
calibration uncertainties. In this case, the oxygen abundance is forced to be lower (higher) 
than the true value with little allowance, because the overall model normalization is tightly 
constrained by the entire continuum and hence cannot vary freely. In fact, the oxygen abun- 
dance and the continuum normalization exhibit strong and negative coupling. In order to 
avoid this problem, we estimated the oxygen abundance by fitting the spectra over a limited 
energy range of 0.5 — 0.8 keV, with all parameters except the oxygen abundance and the oxy- 
gen and continuum normalizations fixed to the best-fit values obtained through the full-band 
fitting. Although this would not be an orthodox method to evaluate a fitting parameter, it 
generally gives a more conservative and unbiased oxygen abundance. 

We firstly fitted the spectra of the three EPIC detectors separately with a single tem- 
perature vAPEC model, in order to study systematic differences among the three EPIC 
detectors. The temperatures derived with PN, MOSl, and M0S2 are kT = 1.26^0:0} keV, 
1.28+[J:|]J keV, and 1.28;[|;|J} keV, respectively. The iron abundances from PN, MOSl, and 
M0S2 are Zpe = 0.23l|]:[5} solar, 0.26;[J:^^ solar, and 0.25^|];!J^ solar, respectively. Thus, the 
three EPIC detectors give consistent results within their statistical errors. In the following 
analysis, we therefore fit the specta of the three detectors simultaneously in order to improve 
the statistics. 

Figure shows the results of the simultaneous fit using the vAPEC and vMEKAL mod- 
els, and Table [T] summarizes the best-fit parameters. Although the obtained parameters are 
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consistent with those derived with ASCA (Paper I), neither one-component model (vMEKAL 
or vAPEC) successfully reproduced the spectra. Since the fit residuals are seen around Fe- 
L line regions, a multi-temperature condition is suggested. Accordingly, we introduced a 
two-component model, consisting of two vAPEC (or vMEKAL) components with different 
free temperatures and normalizations. The two components are assumed to share the same 
metal abundances, and are subjected to the same Galactic absorption. Then, the fit has 
been improved and became acceptable (Figure 5c, 5d, and Table 1), with a < 1% probability 
for this improvement to be accidental. We hereafter adopt the vAPEC code, because it can 
generally reproduce the Fe-L lines better (e.g. Matsushita et al. 2002); our results remain 
essentially the same even if we use the vMEKAL model. With the ASCA data, we derived 
the metal abundance to be 0.3 ± 0.1 using the MEKAL model (Paper I). This is consistent 
with the iron and silicon abundances in Table [T], although the ASCA results are subject to 
larger errors. 

Although we constrained the two vAPEC components to have the same metal abun- 
dances, this may not be necessarily warranted. Therefore, we tentatively relaxed this 
constraint, and found that the iron abundances of the hot and cool components became 
'Z^Fe,High = 0.34 ± 0.02 and 2'jfc,Low = 0.43 ± 0.05, respectively. Therefore, the tied value of 
Zpe = 0.34 ± 0.01 is considered to be appropriate for the hot component, while it may devi- 
ate from the true iron abundance of the cool component by —20%. Although this difference 
exceeds statistical errors, the fit is not improved significantly, because the F-test indicates a 
88% probability of the difference by chance. Furthermore, the oxygen abundances became 
poorly constrained as 2'o,High = 0.18 ± 0.07 and ^o,Low = 0.71 ± 0.68. Considering these, 
we retain our original assumption that the two vAPEC components share the same metal 
abundances. Even when taking into account the possible ~ 20% systematic error in the 
metal abundances of the cool component, our comparison between the metal distributions 
with that of the galaxy light, described in §5, remains essentially unaffected. 

As a next step, we studied spatial variations of the temperature and metal abundances. 
For this purpose, we divided the image into eight concentric annular regions, centered on 
the X-ray emission peak (and hence on the central galaxy). The regions were chosen so that 
each of them contains more than 5,000 photons after the background subtraction. Then, 
the spectra of the eight regions were fitted in the same way as for the spatially integrated 
spectra. Figure [6] shows radial profiles of the temperature, abundances and the reduced 
chi-squared, derived by fitting the individual annular spectra (simultaneous among the three 
detectors) with one-temperature and two-temperature vAPEC models. The temperature 
determined with the one-component model increases mildly from 1.0 keV at the center to 
1.4 keV at ~ 3' (45 kpc), and then returns to 1.0 keV at ~ 10' (150 kpc). Though somewhat 
depending on the modeling (single vs. double temperatures), the oxygen abundance is low 
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and spatially rather constant, while those of the other metals clearly increase toward the 
center by a factor of 3 to 5. The iron and silicon profiles are rather similar, regardless of 
the temperature modeling, although the double-temperature fit yields systematically higher 
abundances than the single-temperature one. In all regions except for the outermost one, 
the two temperature model gives a better reduced chi-squared, and the F-test indicates that 
the probability of these improvements caused by chance is less than 1% in six regions except 
the 5/0 — 7/0 and 10/0 — 14/0 ones. The two temperatures (1.5 and 0.95 keV; Table [1]) 
found with the region-integrated spectra can be understood reasonably as weighted means 
of the two temperatures specified by the eight annular spectra. 

Although the EPIC data thus prefers the two-temperature modeling to the single-phase 
views, this could simply be due to projection effects. We accordingly deprojected the eight 
annular spectra into eight shell spectra using "onion-peeling" method (e.g. Buote et al. 
2003), in which the contribution from the outer shell is subtracted progressively toward 
smaller radii. The emission from the region lying outside the outermost radius (14') was 
estimated by assuming that the spectrum of the outermost annulus extends with its surface 
brightness following the outer jS component derived in §3.1. In the model fitting to the 
deprojected spectra, we set initial parameter values to the best-fit ones obtained with the 
annular spectra, and again fixed or constrained the metal abundances in the same manner as 
the annular analysis. The obtained profiles of the model parameters, presented in Figure U\ 
are similar to those obtained from the annular spectra. The two modelings (single vs. 
two temperatures) are both acceptable essentially in all shells, and their fit goodness no 
longer shows the significant difference. This is due partially to the removal of the projection 
effects, and partially to reduced data statistics. Therefore, the plasma in each shell can 
be represented by a single-temperature coUisionally-ionized plasma model, although local 
two-temperature (or multi-temperature) conditions are not necessarily ruled out. 



4. Near Infra-Red Data Analysis 

In order to estimate the galaxy light profile around NGC 1500, we utilezed the 2MASS 
(Two Micron All Sky Survey) data, which are accessible via NASA/IPAC Exgragalactic 
Database. In the 2MASS catalog, 33 galaxies (with i^-band magnitude less than 14) have 
been found to be lying within 34' of the NGC 1550 nucleus; at the distance of NGC 1550, 
this projected radius corresponds to 500 kpc. The sky positions of the galaxies in 30' x 30' 
region around NGC 1550 are plotted in Figure [3] (right). Among the 33 galaxies, 7 galaxies 
listed in Table [2] are no fainter in i^-band than NGC 1550 by 2.5 magnitude. These galaxies, 
except MCG+00- 11-050 (of which the radial velocity is unknown), have radial velocities 
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within ±500 km s ^ of that of NGC 1550, and are hkely to be members of the "NGC 1550 
group". Three galaxies (IC 0366, UGC 03011, and MCG+00-11-050) of the 7 are within the 
X-ray signal-detected region (14' in radius). 

Figure [8] (red diamonds) shows the 2-dimensional integrated i^-band profile around 
NGC 1550, constructed from NGC 1550 itself and the 33 galaxies. In this figure, the surface 
brightness profile of NGC 1550 is taken into account by using a public 2MASS background- 
subtracted i^-band image of NGC 1550, while the other galaxies were treated as point 
sources. Since this profile must be contributed by background objects, we estimated the 
number of background galaxies expected in the 34'-radius region, based on the number counts 
of galaxies with K < 14 from the entire 2MASS sample (Bell et al. 2003). The /f-band 
magnitude was converted to luminosity assuming the absolute K-hand solar magnitude of 
3.39 (Kochanek et al. 2001) and the same distance to all the galaxies as that to NGC 1550. 
Figure [8] (blue boxes) is the estimated two-dimensionally integrated background i^-band 
luminosity profile. Even at the angular distance of 34', the background is about 5 times 
lower than the integrated i^'-band luminosity of the NGC 1550 system. Within 34' (510 kpc) 
and 14' (210 kpc), the central galaxy carries 39% and 70% of the background-subtracted 
i^-band luminosity, respectively. 



5. Spatial Distributions of Mass Components 
5.1. Mass profiles 

From the parameters of the X-ray emitting gas obtained in the spatial and spectral 
analyses (§3.1 and §3.2), we derived radially integrated mass profiles of three components 
constituting the NGC 1550 system; the ICM mass, the total mass, and the gaseous metal 
mass; the results are shown in Figure [H^. Below, we describe how these results have been 
obtained. 

As a step to calculate the ICM mass profile, we started from the best-fit double-/? model 
of equation [21 We analytically deprojected it into a 3-dimensional emissivity profile as 

where r is 3-dimensional radius from the emission centroid, and {i = 1, 2) are constants 
which are determined uniquely by Sj, r^, and jSi of equation ([2]). Then, the emissivity is 
rewritten as e = A(T, Z)n'^, where A(T, Z) is the cooling function specified by the measured 
temperature and metallicity, and Uc is the electron density. Using this rie, the ICM mass 



- 11 - 



density profile is written as 



A{T,Z) 



(4) 



Here /i = 1.2 is the mean molecular weight (i.e., the mean number of nucleon per electron). 
Values of A(T, Z) were estimated for each shell region, utilizing T and Z measured in § 3.2, 
and referring to tables in Sutherland & Dopita (1993). Finally, we obtained the radially 
integrated ICM mass of Figure [9^, in the standard way, as 



In Figure [9^, the gas mass profile outside 14' was calculated assuming that the emissivity 
profile follows equation with the same temperature and metal abundance as those in the 
10'-14' region. In paper I, the ICM mass inside r = 220 kpc was measured to be 3.4 x lO^^Mo; 
at the same radius, we now obtain a value as Mgas = 7.6 x IO^^Mq. The difference is due 
to the inner j3 component, which is apparent in Figure [Hand included in the XMM-Newton 
analysis. The ICM mass derived here is consistent with that from Chandra (Sun et al. 2003). 

The metal mass density profiles in the ICM were derived by multiplying Pgas('") of 
equation (j4]) with the measured metal abundances. For example, the radial profiles of the 
iron mass density, ppdr), can be derived by multiplying the ICM mass profile with the iron 
abundance profile as 



where Ap^ is the averaged atomic weight of Fe, (ripe/'^H)© = 4.68 x 10~^ is the number 
ratio of iron to hydrogen under one solar abundance (Anders & Grevesse 1989), Zpf,{r) is 
the measured iron abundance, and /h = 0.75 is the hydrogen mass fraction in the ICM. If 
another definition of solar abundance such as Grevesse & Sauval (1998) is adopted, Zpe{r) 
can change, but PFc{f) remains the same, because the change in Zpdr) is canceled out by 
that in ; that is, pFc{r) is independent of the definition of the "solar abundance". 

As to Zpcir) in equation ([6]), we used the results from the single temperature fits to the 
deprojected shell spectra. However, as can be inferred from Figure Ub, the results do not 
change by more than 60% even when the two-temperature fits are employed. 

At r = 210 kpc, the integrated masses of iron, silicon, and oxygen imply 0. 18^005, 
0.271 

ai3' 0-29^a25 solar abundances, respectively, when divided by the integrated gas 
mass. In Figure [9^, iron and silicon exhibit very similar integrated mass profiles, reflecting 




(5) 




(6) 
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similar abundance profiles (Figure [Tt). In contrast, the integrated oxygen mass increases 
somewhat more steeply in Figure [9^ than those of iron and silicon. Obviously, this is because 
the oxygen abundance is radially rather constant, unlike iron and silicon which show the 
central abundance enhancements. In a word, the oxygen mass profile is very similar to that 
of the ICM. 

The integrated total mass profile was obtained assuming hydrostatic pressure equilib- 
rium of the ICM in the gravitational potential, using the formula as 

kTr ( (ilnpgas dXwT 



MUr) = -^\-'^^-^\ , (7) 
/iGm y dmr dmry 

where k is the Boltzmann constant and G is the constant of gravity. In this calculation, the 
temperature T refers to the values measured with the deconvolved spectra. We evaluated the 
term for each region by connecting two adjacent temperature points with a power-law, 
and substituting with the power-law slope averaged on both sides of each point. In 
the present case, — 'j^^ is in the range of 1.2—2.1, whereas is typically < 0.1 (0.5 
at maximum). Therefore, the proper inclusion of the ^-^7 term reduces the estimated Mtot 
up to several tens percent, in those regions where T(r) is increasing outward. Outside 14', 
the temperature was assumed to be radially constant and the same as that in the 10'-14' 
region. At r ~ 210 kpc, we thus obtain the total mass as 2.0 x IO'^^Mq, which agrees with 
the value of 2 x lO^^Afo typically found among galaxy groups (Mulchaey et al. 1996). Also, 
the present result agrees with that in Paper I; 1.6 x 10^'^ M0 at r = 210 kpc. The total mass 
profile obtained with Chandra (Sun et al. 2003) is consistent with our profile (Figure [9^) 
within respective errors. Furthermore, the gas-mass fraction at ~ 400 kpc (half the virial 
radius), Mgas/Mtot = 0.06, is consistent with those observed from galaxy groups at a radius 
of overdensity ~ 500 (Vikhlinin et al. 2006). 

Based on ASCA observation, Dupke & Bregman (2005) reported the presence of a ve- 
locity difference by > 1.5 x 10^ km s~^ in the ICM of NGC 1550. If so, equation ([7]), which 
assumes hydrostatic confinement of the ICM, would become no longer valid. However, a 
claimed velocity difference of 2400 ± 1000 km in the Centaurus cluster, by the same 
authors (Dupke & Bregman 2006), was not found significantly by the latest Suzaku mea- 
surements, with an upper limit velocity of 1400 km s~^ at 90% confidence range (Ota et al. 
2007). Therefore, the reality of the reported bulk motion in the ICM of NGC 1550 is open, 
and hence we retain our original assumption of hydrostatic equilibrium. Even if the bulk 
motion is real, the total mass estimate would change only by a factor of ~ 2 (Ota et al. 
2007), which does not affect our main results. 



The present study requires a quantitative comparison of the X-ray derived 3-dimensional 
radial mass profiles against visible-light distribution representing the stellar mass profile. 
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However, the galaxy distribution obtained in §4 is a 2-dimensional quantity expressed on the 
celestial sphere. Therefore, we have to deproject the galaxy profile into three dimensions, 
or alternatively, project the other three-dimensional profiles (the total mass, gas mass, and 
metal masses) onto two dimensions. We have adopted the latter method, because the galaxies 
aroung NGC 1550 are too sparsely distributed to allow reliable deprojection. Figure[9]D shows 
the projected mass profiles in the integrated form, including the galaxy mass profile. In the 
derivation of the galalxy masses, we assumed i^-band stellar mass-to-light ratio (M/L) of 
all the galaxy to be 1.0 for simplicity, altough M/L in ii'-band varies approximately between 
0.5 and 1.0 according to the color of galaxies (Bell et al. 2003). Thus, the galaxy mass 
profile presented here is essentially the same as the i^-band luminosity profile. 

As can be seen in Figure [Hb, the radially-integrated /^-band mass-to-light ratio of the 
system reaches 75 ± 3 Mq/Lq at 210 kpc. Considering a typical color of galaxies, this 
value is consistent with the 5-band mass-to-light ratio of ~ 380 Mq/Lq (r < 210 kpc) 
derived in Paper I. These values are comparable to those of clusters of galaxies. In fact, the 
mass-to-light ratio of galaxy clusters is typically measured to be ~ 250 Mq/Lq in i?-band 
(Bartelmann & White 2002) and ~ 80 Mq/Lq in K-hand (Kochanek et al. 2003). 



5.2. Iron-mass-to-light ratio profile 

As shown by Figure [9]d, the ICM distribution is slightly more extended than that of the 
total mass (i.e. the gas fraction increasing toward the cluster outskirts), in agreement with 
previous reports (e.g. Markevitch & Vikhlinin 1997). Furthermore, the figure clearly reveals 
that the integrated stellar luminosity profile is significantly fiatter than the other profiles: 
the stellar component is much more concentrated than the other components. The difference 
between the stellar luminosity profile and those of the metal masses is particularly interesting, 
because the metals should have originated from stars. In order to see the difference more 
clearly, we divided the integrated iron mass profile by integrated i^-band luminosity profile, 
to obtain so-called iron-mass-to-light ratio (IMLR) profile (Finoguenov 2000). The radially- 
integrated IMLR profile of NGC 1550, shown in Figure [TOk . drops toward the center by 
about two orders of magnitudes, while fiattens outside ~ 100 kpc. This behavior is in 
agreement with the ASCA result on the Centaurus cluster and Abell 1060 (Makishima et al. 
2001), suggesting that the IMLR decreases towards the center are common among regular 
clusters regardless of their richness. The silicon-mass-to-light ratio (SiMLR) and oxygen- 
mass-to-light ratio (OMLR) profiles also decrease toward the center to the same degree as 
the IMLR. 
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5.3. Density -density plot 

The IMLR and similar quantities compare two relevant mass components in the form of 
ratios between their radially-integrated profiles. This approach is useful to grasp large-scale 
behavior of these components, but may not be suitable to investigate more localized features. 
For this purpose, we devised a new plot which we named density-density plot (DDP). A DDP 
plots densities of a component in individual shells against those of another component in the 
corresponding shells. Although the spatial coordinate is lost in a DDP, we can directly see 
the relation between any two mass components in their differential forms. 

Figure [TUb shows DDPs between the total mass and iron mass (red diamonds), and 
between the gas mass and iron mass (blue squares). The latter DDP indicates that the 
iron mass density increases (toward the center) somewhat more steeply than the ICM mass 
density, just reflecting the iron abundance increase in Figure [7t from ~ 0.15 solar at r ~ 200 
kpc to ~ 0.5 solar at the center. 

Interestingly, we instead find in Figure [TOb a very good proportionality between the 
iron mass density and the total gravitating mass density (red diamonds) except at the very 
center. In fact, the least-square fit with a power law function yields the best relation between 
these two quantities, excluding the innermost data point, as 

ppe = (3.2 ± 1.6) X 10-^ plf^'-'\ (8) 

while the relation between the iron mass density and the ICM mass density becomes 

PFe = (4.7±3.3)xlO-^pic'™- (9) 

In other words, the iron in the ICM traces the total mass distribution (except the innermost 
region) rather than that of the ICM density. In the innermost region (r < 7.5 kpc), the iron 
mass fraction decreases significantly relative to the total mass. This is likely to be related to 
the dominance of the stellar mass of NGC 1550 in this region, although it is not immediately 
clear, at this stage, why the iron mass increases less prominently than the stellar mass. 

Incidentally, equation and equation Q can be combined into a relation between 
the ICM mass density and the total mass density, as picM fx PtoF^^'^^. This is reasonable, 
because the index of 0.77 is close to the typical value of /3 ~ 0.7 (Mohr et al. 1999) found 
among clusters, where P is the same index as employed in equation ([2]). 

In order to investigate the above issue at the innermost region, let us compare densities 
between the galaxy i^-band luminosity and the iron mass. For the same reason as described 
in §5.1, we use the 2-dimensional form, which we named Surface Density-Surface Density Plot 
(SDSDP): a SDSDP plots the surface density of a component in the individual 2-dimensional 
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annular regions against those of another component. Figure [TOb shows a SDSDP between the 
galaxy A'-band luminosity (Sgai) and the iron mass (Spc)- Thus, the innermost two regions 
having > 10®Lq kpc~^ are more optically luminous than the other regions by ~ 2 orders of 
magnitude. This demonstrates the characteristic of the NGC 1550 system as an IXEG. In 
contrast, the iron surface density does not exhibit such a steep increase to the center, so that 
Epe is not proportional to Sgai particularly toward the center. This is another expression of 
the IMLR profile (Figure [TOk): in the central region of the NGC 1550 system, iron in the 
ICM becomes much depleted relative to the stellar light, even though it moderately increases 
relative to the ICM mass density. 

As can be readily inferred, the relation between silicon and the stellar mass is essentially 
the same as that between iron and the stars. Oxygen (though not shown) is even more 
depleted in the central region than iron and silicon, because it lacks the central abundance 
increase. 



Using the XMM-Newton EPIC and 2MASS /T-band data, we derived the radial profiles 
of components constituting the NGC 1550 system; the gas mass, total mass, metal mass, and 
the galaxy luminosity. As seen in Figure [Tt and Figure the metals (oxygen, silicon, and 
iron) are extended to ~ 200 kpc from the center, wherein ~ 70% of the i^'-band luminosity is 
carried by NGC 1550 itself. Most importantly, we have found that iron (as well as silicon) in 
the ICM traces very well the total gravitating mass (eq. [7j), whereas the stellar component 
is significantly more concentrated to r < 15 kpc. This is a non-trivial result, since we would 
naturally expect the metals to trace the stellar mass unless at least either of them evolved 
significantly over the Hubble time in the radial distribution. 

The above finding allows three alternative scenarios; one is that the iron within ~ 100 
kpc has been transported outwards; another is that iron in the central region is contained 
in the central galaxy and hidden from our X-ray view. The other is that the galaxies have 
gradually concentrated toward the center (particularly onto NGC 1550) while supplying 
metals to the surrounding hot gas. 

In the first scenario, metals in the central region must have been transported outward, 
together with the ICM (mostly hydrogen and helium), through diffusion or bulk outflow. As 
discussed in Ezawa et al. (1997; 2001), the diffusion constant D of iron in the hot gas is 
given as 



6. 



Discussion 




(10) 
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Therefore, the iron atoms in the ICM can diffuse at most ~ 10 kpc over the Hubble time; 
this is by far insufficient to explain the present results. In the case of bulk outflow, a large 
amount of hot gas within 100 kpc, comparable in mass to what presently exists there, must 
have been evacuated to realize the IMLR of Figure [TUk and SDSDS of Figure [TUb . Since 
the mass of the hot gas within 100 kpc is 6 x lO^^M© (Figure [9]d), the outflow rate becomes 
~ 50 M0yr~^ if we assume a constant rate over the Hubble time. This is just opposite 
to the case of cooling flows. In other words, the energy needed to drive out the central 
ICM must have been large enough to reverse cooling flows. However, even among medium- 
richness clusters, it is generally considered non-trivial to flnd an energy source which is strong 
enough to suppress cooling flows of this order (e.g. Makishima et al. 2001). Since the depth 
of gravitational potential is comparable to the velocity dispersion of galaxies Ugai, the energy 
required to evacuate the central ICM, -Ecva? becomes enormous, as 

£™~M.„a|., = 3xlO"(^35^)%rg. (11) 

where Mcva ~ 1 x 10^'^ Mq is a typical ICM mass to be evacuated. 

The most popular candidate for driving the outflow is activity of the central galaxy 
nucleus. Actually, NGC 1550 hosts a weak radio source of which the flux density is 16.6 ± 1.6 
mJy at 1.4 GHz (Condon et al. 1998). However, for the estimated kinetic energy (power 
times age) provided by an AGN to become comparable to equation ( |TT1) . the object needs 
to be as strong a radio source as Cygnus A (e.g. Kino et al. 2005); the present radio 
luminosity of the NGC 1550 nucleus falls 6 orders of magnitude in short. Of course, we cannot 
exclude an extreme possibility that the nucleus remained highly active (either continuously 
or intermittently) till a near past, and then declined at present. In such a case, we would 
observe large (~ 100 kpc) cavities in the ICM, but there is no such evidence. Alternatively, 
the radiation pressure from the AGN may have caused the outflow. If we assume the ICM 
within 100 kpc from the center has been pushed away, the luminosity of the AGN has to be 

Lagn= (^^j = 4 X 10^6 erg s-\ (12) 

where E^va is the evacuation power averaged over the Hubble time, and tt = 2 x 10"^ is the 
cross section for Thomson scattering. It is unlikely that the NGC 1550 nucleus has been so 
luminous for the Hubble time. Yet another weak point of the scenario is that it needs some 
flne tuning mechanisms to make ppe proportional to ptot as indicated by equation and 
the DDP of Figure [TUb . Thus, the bulk outflow scenario is also unlikely. 

The second scenario assumes that the IMLR decreases toward the center because a 
signiflcant fraction of iron in the central region has been taken inside the stars in NGC 1550, 
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thus increasing the stellar metallicity of NGC 1550 compared to those of the other peripheral 
galaxies. This "hidden" amount of iron becomes comparable to that in the hot gas within 
100 kpc, that is, 2 x IO^Mq (Figure [Hb)- For comparison, the stellar mass in NGC 1550 is 
estimated to be 2 x IO^^Mq, assuming again the i^-band stellar mass-to-light ratio to be 
1. Therefore, the hidden iron should increase the mean stellar metal abundance by ~ 0.5 
in solar units. In this estimate, we can ignore the amount of iron confined in interstellar 
matter (ISM) in these galaxies, because the mass of ISM is at most 10% of the stellar mass 
(Matsushita 2000). From a color- metallicity relation studied by Barmby et al. (2000), the 
stellar metallicity of NGC 1550 is estimated to be 0.5 — 0.7 solar abundance, using U,B,V,J, 
and K magnitude. On the other hand, Jordan et al. (2004) showed that the metallicity of 
early-type galaxies correlates well with the absolute magnitude My- From this correlation, 
the metallicity of galaxies having the same My as NGC 1550 is estimated as ~ 0.8 solar. 
Since this is not much different from the value estimated from the color of NGC 1550, there is 
no reason to presume that NGC 1550 has locked a much larger fraction of its metal products 
into its stars than other elliptical galaxies with simillar l^-band luminosities. 

The remaining scenario assumes gradual concentration, or infall, of the member galaxies 
toward the system center. This view is suggested by the present and other (e.g., Makishima 
et al. 2001) results, that the stellar mass in groups and clusters are generally much more 
concentrated, not only than the metals in the ICM, but also than the ICM itself and the 
total mass (Figure [9b). The former two scenarios would not easily explain these differences 
among the metals, stars, and the total mass (dominated by the dark matter). To be more 
quantitative, let us assume, as a simplest case, that the system initially consisted mainly of 
~ 10 galaxies each with i^' ~ 11, which is comparable to IC 0366 but an order of magnitude 
less luminous than the present-day NGC 1550. If only one of them was initially at r < 50 kpc 
while the others were outside, the IMLR at this radius was initially an order of magnitude 
higher than it is at present, and hence the IMLR profile at that time was approximately 
flat as would be naturally expected. If these galaxies have all merged into the central single 
galaxy (i.e., the present-day NGC 1550) after their metal supply was mostly over, we can 
explain both its present optical luminosity, and the central decrease in the IMLR profile. In 
other words, NGC 1550 is indeed suggested to be a remnant of such a merger process, as 
implied by the nomenclature of "fossil group" (Jones et al. 2003). 

As to the mechanism which drives the suggested galaxy infall, a promising candidate is 
dynamical pressure exerted onto galaxies when they swim through the ICM. According to 
Sarazin (1988), this effect will cause the kinetic energy of a galaxy, moving with a velocity 
V, to decrease at a rate of 

dK/dt= -Tipg^v^Rl, , (13) 
where i?int is "interaction radius" with which the galaxy interacts with the ICM and displaces 
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it. Then, the galaxy will lose its kinetic energy on a time scale of 

^ K ^ r»gai 
^ " \dK/dt\ ~ 27rpgasyi??nt ' 

where mgai is the galaxy mass. 

As discussed by Makishima et al. (2001), the value of i?int is estimated to be comparable 
to the visual galaxy size, because the in-flowing ICM will be intercepted by magnetic fields 
anchored to the galaxy. Assuming hence -Rint = 10 kpc, and employing plausible values for 
the other parameters, equation (fT4l) yields 

r = u( 1 ( \ ^ ( ) ^ ( Gyr (15) 

Vl X lOiiMo; VlO-3 cm-3; V400 kms-i/ VlOkpcy ^ ^ ^ 

which is of the same order as the Hubble time. Then, the galaxies moving in dense ICM 
will fall toward the center, merge together, and change their morphology (Makishima et al. 
2001), on the Hubble time scale. Importantly, equation ( fTSi) predicts less massive galaxies 
to fall more rapidly (on condition that i?int depends only weakly on mgai), in agreement with 
the idea that small galaxies which may have existed in the system have mostly merged into 
NGC 1550. This mass dependence is opposite to that of the gravitational dynamical friction, 
which works among galaxies moving in vacuum and causes more massive galaxies to sink to 
the potential center. 

The scenario described above is not considered specific to the NGC 1550 system, since 
the dynamical pressure cannot be avoided as long as galaxies are moving through the ICM. 
Employing the empirical relation among the cluster mass M, the X-ray luminosity L^, and 
the temperature T as (x and M oc T^'^ (e.g. Arnaud et al. 2005), The ICM density 
rigas in a cluster is found to scale as oc T'^/^; here we also utilized the virial relation as 
M/R oc T, and the luminosity relation as oc n^^R'^T^^^ , where R is the cluster radius. 
Then, with v oc T^/^, equation (IMI) predicts r oc T~^, so that richer systems are expected 
to have somewhat shorter merging time scales. In fact, the central decrease in the IMLR 
profile has been observed from richer systems as well (Finoguenov et al. 2000; Makishima et 
al. 2001; Kawaharada 2006), including the Centaurus cluster and Abell 1060. Interestingly, 
in Figure 9 of Finoguenov et al. (2000), the IMLR profiles in HCG 51 and HCG 62 are 
flatter than those of the other groups and clusters. This support the galaxy-infall scenario 
indirectly, because a compact group is considered to be a galaxy group in which galaxy 
mergers are proceeding. Furthermore, the ICM dynamical pressure might also explain some 
evolutionary effects observed from rich clusters, such as the long-known Butcher- Oemler 
effect (Butcher & Oemler 1978), and the significant evolution of galaxy morphology in dense 
cluster environments (e.g., van Dokkum et al. 2000; Goto 2005). 



(14) 
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It is known that galaxy groups exhibit rather large scatter in their gas richness (Osmond 
& Ponman 2004), with spiral- rich ones tending to be particularly gas poor (Ota et al 2004). 
From equation (IT^ . we may then speculate that a galaxy group borne as a relatively gas-rich 
system is subject to a fast merger process, while it remains otherwise an ordinary group. If 
this view is correct, an IXEGs could be a system in which baryons remained for some reason 
mostly in the hot gas phase, with a relatively lower galaxy mass. The NGC 1550 system 
agrees with this view, because it is as gas-rich as ordinary galaxy groups in terms of the 
gas-to-total mass ratio (Paper I), while it is optically comparable to a single galaxy. 

We acknowledge support from the grant for young researchers from Japan Society for the 
Promotion of Science (JSPS), and the grant from Special Postdoctoral Researchers Program 
of RIKEN. 
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Table 1. Results of the model fitting to the MOS and PN spectra within 14' of the 

NGC 1550 nucleus. 









abundances (solar) 






Model 


Temperature (keV) 


Ratio'' 


Fe 


Si 


O 


xVdof 




vMEKAL 


1.32 ±0.01 




0.29 ± 0.01 


0.31 ±0.02 


0.13 ± 0.03 


2012/1423 = 


1.41 


vAPEC 


1.27 ±0.01 




0.24 ±0.01 


0.32 ±0.02 


0.13 ± 0.03 


1813/1425 = 


1.27 


2vMEKAL 


1.50 ±0.02 / 0.87 ±0.01 


5.91 


0.42 ± 0.01 


0.45 ±0.02 


0.24 ±0.04 


1640/1421 = 


1.15 


2vAPEC 


1.48 ±0.02 / 0.95 ±0.01 


4.08 


0.34 ±0.01 


0.40 ± 0.02 


0.18 ± 0.03 


1666/1423 = 


1.17 



^Ratio of the hot-component normalization to that of the cool component 



Table 2: Major candidate of member galaxies in the NGC 1550 system. 



Name 


(arcmin) 


K (mag) 


Vr (km s ^) 


NGC 1550 





8.774 


3714 


IC 0366 


3.1 


11.104 


3736 


UGC 03011 


12.0 


10.891 


4207 


MCG+00-11-050 


12.7 


10.840 




UGC 03008 


17.0 


9.687 


3251 


CGCG 393-007 


17.1 


10.308 


3912 


CGCG 393-008 


23.6 


9.955 


3973 


UGC 03006 


33.2 


9.581 


3664 
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Fig. 1. — MOSl 10 — 12 keV count-rate histograms of the NGC 1550 observation (dark blue) 
and the template background (orange) with 100 sec integration, both normalized to their 
exposures. 
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Fig. 2.— 0.5-9.0 keV Raw spectra of PN (black), MOSl (red), and M0S2 (green) extracted 
from the 0' — 14' region. For comparison, estimated background spectra of PN (blue), MOSl 
(cyan), and M0S2 (magenta) are also plotted. 
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Fig. 3.— (left) The 0.5 - 4.5 keV EPIC MOSl surface brightness (contours) of NGC 1550, 
including background, overlaid on an optical image of Digitazed Sky Survey (greyscale) 
plotted on J2000 {a, 6) coordinates. The X-ray image is smoothed by a gaussian filter 
with a = 20", but not corrected for vignetting. The contours are logarithmically spaced, 
with a step of factor 1.4. (right) The same optical image of Digitazed Sky Survey as left 
panel, where also plotted are positions of galaxies around NGC 1550 identified in the 2MASS 
catalog. Colors specify their i^-band magnitudes as 8 < -R' < 9 in red, 9 < < 10 in green, 
10 < -fC < 11 in blue, 11 < < 12 in magenta, and 13 < K in cyan. Galaxies with 
12 < < 13 are not present in this sky region. 
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Fig. 4. — The 0.5 — 4.5 keV background-subtracted and vignetting-corrected X-ray surface 
brightness of NGC 1550 created with the MOSl data (black crosses), shown together with 
the best-fit double-/? model (red). The data were azimuthally averaged around the NGC 1550 
nucleus. The constituent two /3 components are shown in blue and green. 
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Fig. 5. — Response-inclusive 0.5—5.0 keV PN (black), MOSl (red) and M0S2 (green) spectra 
from the 0' — 14' region, simultaneously fitted with; (a) an absorbed vMEKAL model, (b) an 
absorbed vAPEC model, (c) a two-temperature vMEKAL model, and (d) a two-temperature 
vAPEC model. In the two temperature fits, the metal abundances are tied between the two 
components, while the temperatures and normalizations are left free. 
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(b) NGC 1550 Normalization 
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d) NGC 1550 Oxygen Abundance 
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(e) NGC 1550 Reduced chi squared 
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Fig. 6. — Radial profiles of the X-ray spectroscopic parameters, derived with 8 annular 
spectra, (a) The temperature derived with the single (black) and double (red and blue) 
vAPEC fits, (b) The vAPEC normalizations (proportional to emission measure), with the 
same color definitions as in panel (a), (c) The iron (magenta) and silicon (green) abundances 
obtained from the single (solid diamonds) and double (dashed diamonds) vAPEC fits, (d) 
The oxygen abundances, with the same definitions of solid and dashed diamonds as (c). (e) 
The reduced chi-squared from the single (black) and double (red) vAPEC fits. 
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(I)) NGC 1550 Normalization 
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(d) NGC 1550 Oxygen Abundance 
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Fig. 7.- 



The same as Figure [6l but for the deprojected spectra. 
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Fig. 8. — A two dimensionally integrated i^-band luminosity profile around NGC 1550 (red 
diamonds), compared with that of the estimated field galaxies (blue squares). The latter is 
calculated assuming that all the galaxies are at the same distance as NGC 1550. 
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(a) NGC1550 3D mass profile 
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(b) NGC1550 ED mass profile 
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Fig. 9. — Integrated radial profiles of various mass components in the NGC 1550 system, 
shown in 3-dimensional forms (panel a) and projected forms (panel b). Black crosses indicate 
the total mass calculated by equation ([7j), while blue crosses the total gas mass obtained 
by equation ([5]). Small cyan diamonds are the galaxy mass derived in §4, shown only in 
the projected form. Metals in the gas are specified by colors; O (red diamonds). Si (green 
squares), and Fe (magenta circles). 
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(b) N1550 DDP 
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Fig. 10. — (a) Radially-integrated i^'-band iron-mass-to-light (magenta), silicon-mass-to- 
light (green), and oxygen-mass-to-light (red) ratio profiles of NGC 1550. For convenience, 
the angular separation is converted to physical distance, assuming that all the galaxies are 
at the same distance as NGC 1550. (b) Density-density plots between the total mass and 
iron mass (red diamonds), and between the gas mass and iron mass (blue squares). Two 
dashed lines indicate 0.1 and 1.0 solar abundances of the iron mass. Green dotted line 
shows the best fit power-law relation between the iron density and the total mass density, 
Pfc = 3.191 X 10~^ Ptot''^) excluding the innermost data point, (c) A surface density-surface 
density plot between the galaxy /T-band luminosity and iron mass (red diamonds). The 
dashed line represents oc Sgai for reference. 



